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ABSTRACT 

The nature of the supernova leading to the Crab Nebula has long been controversial because of 
the low energy that is present in the observed nebula. One possibility is that there is significant 
energy in extended fast material around the Crab but searches for such material have not led to 
detections. An electron capture supernova model can plausibly account for the low energy and the 
observed abundances in the Crab. Here, we examine the evolution of the Crab pulsar wind nebula 
inside a freely expanding supernova and find that the observed properties are most consistent with a 
low energy event. Both the velocity and radius of the shell material, and the amount of gas swept up 
by the pulsar wind point to a low explosion energy 10 50 ergs). We do not favor a model in which 
circumstellar interaction powers the supernova luminosity near maximum light because the required 
mass would limit the freely expanding ejecta. 

Subject headings: ISM: individual objects (Crab Nebula) — supernovae: general — supernovae: indi¬ 
vidual (SN 1054) 


1. INTRODUCTION 

The identification of the supernova type of SN 1054, 
the event leading to the Crab Nebula, has been an en¬ 
during mystery. The amoun t of m ass in the Crab fila¬ 
ments, 4.6T1.8 M© ( Fesen et al.|1997l), and their typical 
velocity, 1260 km s -i (Temim et al. }2006 ), imply a ki¬ 
netic energy of ~ 7 x 10®" 

ergs, considerably less than 
the ^ 1 x 10 51 ergs inferred for typical core collapse su¬ 
pernovae. One possible reason for the discrepancy is that 
what we view as the Crab Nebula is only the inner region 
of a more extended supernova event, which has a normal 


energy, ^ 10 51 ergs (Chevalier 1977). A constant den¬ 


sity freely expanding region with this energy and mass of 


6 M© would extend to ~ 5300 km s i . However, there 
is no sign of such a fast ejecta envel ope int eracting wi th 


the surrounding medi um at X-ray (Seward et al. 2006) 


or radio wavelengths (Frail et al. 1995| ), which suggests 


that the density of the surrounding medium is low and 
the fast ejecta would be in free expansion. In that case, 
the freely expanding ejecta are heated and p hotoionized 
by the radiation from th e Crab Nebula (Lundqvist et al. 
1986; Wang et al.||2013 ), but they have not been clea rly 


detected in either emission lines (Fesen et al. 
in absorption to hig h velocities alo ng the line o: 
the central pulsar ISollerman et ah[ 2000 


Tziamtzis||2012 Lundqvist et al.||2015fr 


1997) or 
sight to 


Lundqvist & 


If the surrounding ^ 10 5i erg supernova remnant is 
absent, SN 1054 must have been a low energy and rel¬ 
atively low mass event, and there is the oretical support 


for such a view. Nomoto et al. (1982) found that the 


abundances in the Crab INebula pointed to a progenitor 
mass in the range 8 — 10 M© and a probable electron 
capture supernova. Models of electron capture super¬ 
novae have shown that neutrino driven explosions occur, 
but with a relatively small energy > 1 x 10 50 ergs (Ki- 
taura et al. 2006). On the observational side, low energy 
(and low velocity) super novae have been discovered, in¬ 
cluding IIP events such as SN 1997D, with an estimated 


energy of 10 50 ergs (Chugai & Utrobin||2000). However, 
SN 1997D had a peak absolute magnitude of —14, con¬ 
s idera bly faint er than SN 1054 at maximum. Chugai fe 
Utrobin (2000) suggested that SN 1054 was a low energy 
supernova, with the lumin osity provided by circumstellar 
interaction 

994W and SN 201 lht are likely exam¬ 


ples of su ch events and are poten tial analogs to SN 1054. 
However, Tominaga et al. (2013) followed the presuper¬ 
nova evolution as well as the explosion of an electron 
capture supernova, finding that a peak luminosity com¬ 
parable to SN 1054 could be produced without strong 
circumstellar interaction; the high peak luminosity is re¬ 
lated to the large extent of the progenitor star. 

In the circumstellar interaction scenario for the peak 
emission from SN 1054, the dense filaments in the Crab 
contain circumstellar material and the medium around 
the observed nebulae is the initial surroundings, not su¬ 
pernova ejecta. An alternate view, which we investigate 
here, is that the nebula is surrounded^by freel y expand¬ 
ing ejecta, as in the model of Chevalier (1977), but we 
allow for a low energy supernova. Basic constraints on 
the model are discussed in Section 2. Pulsar wind bubble 
models are presented in Section 3 and the conclusions are 
in Section 4. The observations of the Crab are spread 
over the past decades. Here we take 1994 as a typical 
time of observation so that the age of the Crab Nebula 
is 9 40 yr. The distance to the Crab is taken to be 2 kpc 
(Trimble| |1973 ). 


2. OBSERVATIONAL CONSTRAINTS 

The view taken here is that the outer boundary of the 
observed Crab is the interface of the pulsar wind neb¬ 
ula with freely expanding ejecta to the outside, as in the 
model of Chevalier (1977). The fact that the synchrotron 
emitting nebula fills the region of thermal emission is 
highly suggestive of the outer boundary being the bound¬ 
ary with freel y expanding ejecta. Additionally, San krit fc 
Hester] (]1997h found that the emission of C IV and Ne V 
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Fig. 1.— Schematic diagram of the model for the Crab Nebula 
discussed in the text. 


at the outer edge of the Crab is consistent with shock 
emission where the pulsar bubble is driving into freely 
expanding ejecta, but is inconsistent with a photoioniza- 


to come from a “skin” to the Crab (|Hester 

to 

o 

o 

00 

|. The 

basic model presented by bankrit fc Hester 

(T997) has 


a shock velocity of v s h = 150 km s -i and a preshock 
density of po = 12.7 mn cm -3 , where m# is the proton 
mass. The model requires that there be a cooling shock 
wave to produce the line emission, although there may 
not be such a shock front around the whole nebula (Loll] 
et al.||2013 ). The situation is illustrated in Figure fl] 


The outer layers of the main body of the Crab Nebula 
are expanding at 1700 — 1800 km s -1 . If this is cooled 
postshock gas, the free expansion velocity at the outer 
edge is 1550 — 1650 km s -1 if there is a 150 km s -1 
shock front. There is evidence for freely expanding gas 
extending to ~ 2500 km s -1 , well beyond the velocity 
of the Crab shell material. One piece is the spectrum of 
the central pulsar, whi ch shows C IV absorpt i on extend¬ 
ing to ~ 25 00 km s -1 ( Sollerman et al.||2000 Lundqvist 


et al. 2015). The absorption does not show the broad 
hat trough that would be expected if the Crab were sur¬ 
rounded by the extended envelope of a normal 10 51 erg 
supernova. A pulsar wind nebula is expected to have 
higher velocities in the shocked wind close to the pulsar, 
but there has been no detection of high velocity thermal 
gas in this region; presumably the gas has been swept 
out by the flow of relativistic fluid. 

The other piece of evidenc e comes fr om the “jet” on 
the north side of the Crab. Rudie et al. (2008) found 
that the emission along the sides of the jet is consistent 
with freely expanding gas that has been moving with 
approximately constant velocity since 1054. The outer 
tip of the jet extends to ~ 2500 km s -1 (Rudie et al. 


2008 Black fc Fesen|2015 ). A plausible origin for the jet 


is that relativistic fluid from the main body of the Crab 
has pushed out and moved aside freely expanding gas 
outside the Crab. In this way, o utwa rd acc elerati on of 
the gas could be relatively small. Rudie et al. ( 2008|) refer 
to this picture as the aneurism model. |banknt fc Hester 


(1997) advocate such a model to explain the [O IIIJ skin 


at the sides of the jet. These considerations suggest a 
model in which the density is declining outside the main 
body of the Crab and freely expanding gas extends to 
> 2500 km s- 1 . 

Another constraint comes from the acceleration of the 
optical filaments determined by proper motions. The an¬ 
gular distance from the center of expansion divided by 
the proper motion gives the year of explosion on the as ¬ 
sumption of constant velocity: 1140 ± 15 QTrimble|T 968), 
11201 =b 7 (j Wyckoff fc Murray||1977|) , and 1130 ± 16 (|Nu- 
gent 1998), corresponding to Al' of 86 db 15, 66 ± 7, and 
76 ±16 yr, respectively, where AT is the difference be¬ 
tween the actual explosion date (1054) and the constant 
velocity explosion date. These results are for optical fil¬ 
aments that are spread around the extent of the Crab 
Nebula. A complication with the Crab is that the main 
shell extends over a range of velocities. The o uter edge 
of the nebula is moving at 1700 km s 1 ( |Cla rk eb 
ah]|1983 Sollerman et al.||~2000 ), but the b rightest emis- 
sioryshows shell expansion at 1260 km s -1 (Temim et al. 
2006) and many of the observed knots are probably in 
this structure. The free e xpansion velocity at the outer 
edge is ~ 1600 km s -1 . (Bietenholz et al.| (|l99l|) have 


examined the expansion of the outer radio synchrotron 
emission, obtaining an explosion date of 1245 ±92, or 
AT = 191 ± 92 yr. In the picture of filament formation 
by Rayleigh-Taylor instability, there is the expectation 
of higher acceleration in the outer parts of the nebula, 
and there is some evidence for such higher acceleration. 

For the model in which the pulsar nebula is expanding 
into freely expanding ejecta, the shock velocity can be 
expressed as v s h = V s — ( R/t a ), where V s = dR/dt is 
the velocity of matter at the outer edge of the nebula in 
the rest frame and t a is the true age of the Crab. For 
material at the outer edge, 


AT = t — — = _—_ 

a V, 1 + (R/taV8hY 


(1) 


For v s h = 150 km s -1 , R/t a = 1600 km s -1 , and t a = 940 
yr, we have AT = 81 yr. The acceleration needed to 
produce the shock velocity in the freely expanding ejecta 
is roughly consistent with the observed acceleration. 

The SN 1054 progenitor is estimated to have an initial 
mass of 8 — 10 based on the observed abundances 
(Nomoto et al. 1982). Allowing for 1.4 M 0 going into 
the central neutron star and 1 M 0 , or more, going into 
presupernova mass loss, we have an ejecta mass, Msn , of 
6.6±1 M 0 , or less if there is more mass loss. The current 
Cra b Nebula has a mass M sw ~ 4.6 ±1.8 M 0 (Fesen et 
aL[l997). The ratio of swept up mass to total supernova 
mass is then M sw /Msn ~ 0.70 ± 0.30. If the observed 
nebular mass is composed of supernova ejecta swept up 
by the pulsar wind nebula, as in the model advocated 
here, much of the ejecta are likely to have been swept 
up. 


3. EVOLUTION OF THE CRAB NEBULA 

There are several components of the evolutionary 
model. In the presupernova stage, there is mass loss 
that provides the outer target for the expanding super¬ 
nova. Within days after the supernova, the supernova 
ejecta approach free expansion. The outer parts interact 
with the presupernova mass loss, while the inner parts 
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Fig. 2. — Density distribution of the freely expanding supernova 
gas for an approximate power law structure and for a more de¬ 
tailed model of an exploded red supergiant star. The indices S 
and p represent the small velocity and large velocity asymptotes, 
respectively, for the detailed model. The approximate model with 
power law index m is that given by equation (pi. The assumed 
ejecta mass is 5 Mq, the supernova energy is 10°° ergs, and the 
age is 940 yr. 

interact with an expanding pulsar wind bubble. 


3.1. Approximate Model 

Our aim is to examine the evolution of the Crab Neb¬ 
ula, considering the dependence on the supernova energy. 
We begin with the supernova explosion, which results in 
an inner flat and outer steep density profile. The ejecta 
are expected to reach free expansion within a few days 
and we approximate the profile by 

( At~ 3 for V <V L 

P= \At- 3 ^y 7 iovV>V t , (2) 


where A is a constant, V = r/t is the velocity of freely 
expanding gas, and V t is the transition velocity from a 
flat to a steep density profile. A comparison to the more 
accurate profile used in Section 3.3 is shown in Figure 
[2| The mass and energy of the supernova are given by 
AfgN = 7ttAV^/3 and Esn = 77rAv^ 5 /5, so that V t = 

(SE’sn/SMsn ) 1 / 2 = 1.29 x 10 3 E 1 5 ( ) 2 M~ 1/2 km s" 1 , where 
E 5 o = -EsnAIO 50 ergs) and M 5 = M SN /5 M©. Also, 


A = 6.34 x 10 8 E~ 0 3/2 M 3 /2 g s 3 cm" 3 . 

On a timescale of 100 ’s of years, the pulsar bubble 
expands into the freely expanding ejecta on the inside. 
For the Crab pulsar, we take a period P = 0.0331 s, 
P = 4.2 x I Q " 13 s s -1 , and a constant braking index 
n = 2.5 (Lyne et al. 1988). Using 


P = Po 



( 3 ) 


where Po is the initial period and r is a constant 
timescale, and the above parameters, we find r = 725 
yr for t = 940 yr. With a neutron star moment of 
inertia of 10 45 g cm 2 , the current spindown power is 
E = 4.6 x 10 38 erg s -1 and it evolves as 

E = E 0 (l + Aj " _1 , (4) 


where Eg is the initial power. Assuming the observations 


were made at an age of 940 years, we have Eo = 3.2 x 
10 39 erg s -1 . The average power input up to time t is 


En = 


(n - l)E 0 r 
2 1 


1 - I 1 + - 


-2/(n— 1)' 


( 5 ) 


For the Crab pulsar parameters, E a = 0.39Eq = 1.2 x 
10 39 erg s -1 at the present time. 

The expansion of the pulsar bubble is determined by 
the equ ations o f conser vation of mass, momentum, and 
energy (Chevalier 1977): 

dMsw = 4t tR 2 P (V s - V), ( 6 ) 


M, 


dt 

■w-<•*[* 

d(A-KR A Pi] 

dt 


P(V S - V) 2 ] 


= ER , 


( 7 ) 


(8) 


where Pi is the interior bubble pressure and it has been 
assumed that the pulsar bubble expands adiabatically 
with an adiabatic index 7 = 4/3. The thin shell approx¬ 
imation is made in these equations. 

The pulsar bubble initially expands in the constant 
density inner region of the ejecta. If the pulsar power 
can be regarded as consta nt, the expansio n is self-similar 
and the bubble radius is (Chevalier 1977) 


R = 


875 


396 


V[Ea\ 


1/5 


Msn J 


t 6 / 5 . 


(9) 


This evolution implies that the time to sweep out the 
ejecta to velocity V t is 


tt 


132 Esn 

1751 ^ 7 ' 


( 10 ) 


Taking E a for the Crab pulsar, we have t t = 1940^50 
yr, so the shell is still primarily inside of the density 
inflection point. Equation § can be used to find the 
velocity of the shell at the present age of the Crab 

Vs = Ij= 1.28 x 10 3 EU 10 M~ 1/2 E 1 a il km sT (11) 

The velocity is less than that at the outer edge of the 
Crab Nebula, but is comparable to that deduced from 
observations of bright filaments in the Crab. An impor¬ 
tant point is that a supernova energy of 10 5 ° ergs gives an 
outer velocity in approximate agreement with the Crab 
Nebula dense filaments, while 10 51 ergs gives too high a 
velocity. 

In this model, the ratio of swept up to total ejecta mass 
is 

/ • \ 3/5 

= 0.68 ( M J = 0.33E 5 -°- 6 E ° 3 6 9 (12) 


for the age of the Crab Nebula. This value is at the low 
end of what is allowed by estimates of the observed swept 
up mass and the mass of the progenitor star. Going to 
a supernova energy of 1 x 10 51 ergs would decrease the 
predicted ratio by a factor of 4 to well below the observed 
value. 
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In this simple model, V s = (6/5 )(R/t a ), so that 

at = t a - ~ ^ = 157 yr, 

V* 6 


(13) 


independent of other parameters. The result is larger 
than that deduced from optical observations. If there are 
the dependences E oc t~ s and central density A oc r -5 , 
dimensional analysis shows that the 1/6 factor in equa¬ 
tion ( p~3] ) becomes (1 — s + 5)/(6 — s). The density is 
expected to drop with radius, giving larger acceleration 
and an increasing discrepancy between the theory and 
observation. The power input E is expected to drop with 
time, decreasing the acceleration and improving the dis¬ 
crepancy between the theory and observation. The more 
detailed model presented in Section 3.3 has the capability 
to give better agreement with the observations. 

The density in the flat part of the supernova 
density profile at an age 940 yr is pf = 2.4 x 

10~ 23 M^ 2 E 3 ^ 2 g cm -3 . If the free expansion velocity 
at the outer edge of the Crab is 1600 km s -1 , as discussed 
above, the density at that point is 0.54 x 10 -23 g cm -3 , 
which is somewhat belo w the preshock density of 2.1 x 
IQ- 23 g cm -3 inferred by Sankrit & Hester (1997) in their 
radiative shock model. An energy of 1.6 x 10 5U ergs would 
improve the agreement, but with an energy of 1 x 10 51 
ergs the highest ejecta density is 0.077 x 10 -23 g cm -3 . 
The radiative shock model is more consistent with a low 
energy explosion. The observed shell is only partial, 


ter 


at lower shell velocities (Sankrit & Hester 1997 Hes 


2008), as would be expected for the density profile 


proposed here. We note that the outer edge is in the 
region of declining density here, whereas it is within this 
region in the discussion of the radial expansion because 
the outer edge of the Crab Nebula is substantially outside 
of the dense filaments. This issue is further considered 
in Section 3.3. 

3.2. Interaction with Dense Mass Loss 


In the SN Iln-P model of Smith (2013), the super¬ 
nova interacts with dense circumstellarnatter, with 
mass Mcs, on a timescale of ^ 10 2 days. The in¬ 
teraction may be responsible for the optical light from 
the supernova near maximum. In the presence of the 
circumstellar material, the density profile in equation 
© has a cut-off, at some velocity V c . All the mate¬ 
rial with speed higher than V c interacts with circum¬ 
stellar material and forms a dense shell with velocity 
V c . We initially assume that the presupernova velocity 
of the initial circumstellar material is negligible, com¬ 
pared to the supernova expansion. Under this assump¬ 
tion, the radial momentum of matter with velocity > V c 
is P r (V > V c ) = [M(V > V c ) + M CS } x V c - Assuming 
V c > V t , we find P r (V > V c ) = M(V > V c ) x 4U c /3, 
leading to M(V > V c ) = 3Mcs after the circumstel¬ 
lar matter has been swept up. The total mass of the 
shell, including circumstellar and ejecta components, is 

M Sh ell=^M C S- 

For V c > Vt we find 

v -(^r v - <i4 » 

so the condition for V c > V t is Msn > 7 Mcs • When 


TABLE 1 

Supernova to circumstellar mass ratio 


Vc/Vt 

M S n/M C s 

0.25 

0.33 

0.5 

0.98 

0.75 

2.5 

1.0 

7.0 

1.5 

35 

2.0 

112 


V c < V t , an integral shows that the relation is 

-1 /TT \ 3 


Mcs 

Msn 


i + I(T 

7 \V t 


(15) 


Values of Msn /Mcs for various values of V c /Vt are given 
in Tabled] 

We have assumed here that the swept up shell of ejecta 
and circumstellar medium is thin because the internal 
energy generated by the interaction is carried away by 
radiation. This is a reasonable assumption in cases where 
the radiation is inferre d to carry a signific ant part of the 
supernova energy (e.g.,|Chugai et al.|2004|). The radiated 
energy is then 


Erad = AE 1 = Esn{V > V c ) ~ -MshellVc , (16) 

where Esn{V > V c ) is the energy in the freely expanding 
ejecta with velocities > V c . If V c > Vt we find 


Erad = k EsN (y t 


-2 


25 E% 


63 M S n 


SN v -2 
V r % 


(17) 


so that E rad = 4.0 x 10 49 E 2 0 M 5 1 V c3 2 ergs, where V C 3 = 
V c /1000 km s _1 . The value of E rad for SN 1054 is not 
well determined from the Chinese observations. The light 
curve is roughly consistent with a normal Type II super 


nov a, so that E rad ~ (1 — 3) x 10 49 ergs. [Mauerhan et 
al. (2013) estimated a radiated energy of (2 — 3) x 10 4y 
ergs for the Type Iln-P SN 201 lht. If V c = we have 
Erad = 2.4 x 10 49 E 5 o ergs, which suggests that V c is close 
to V t or 1.29 x 10 3 km s -1 for the standard parameters. 
However, the profile of freely expanding gas is cut off at 
V c and we argued in Section 2 that there is freely ex¬ 
panding gas over the velocity range 1600 — 2500 km s -1 . 
The high velocity of free expansion suggests that V c is at 
least twice Vt, which drops E rad by a factor of 4. There 
is thus a discrepancy between the Type Iln-P model for 
SN 1054 and the evidence for freely expanding ejecta in 
the Crab Nebula. 


In the model of Tominaga et al. (2013), the progen¬ 
itor is a super asymptotic giant (SAGE!) star with an 
estimated mass loss rate of 10 -4 M e yr -1 and a wind 
velocity of 29 km s -1 . In their discussion of SN 1054 


as an electron capture event, |Moriya et al. (2014J) took 
the SAGB wind velocity to be 10 km s -i . The wind du¬ 
ration is ^ 10 4 yr, so it extends to (3 — 9) x 10 17 cm, 
or (0.1 — 0.3) pc in radius. Unlike the case of a Type 
Iln-P supernova, the forward shock wave is nonradiative 
from an early age, although the reverse shock is radiative 
while the interaction region in is the dense wind, for the 
standard parameters. For the supernova density profile 
given in equation the interaction shell moves through 
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the wind on a timescale of 200 yr. The velocity of freely 
expanding gas at that point is 1460 km s _1 . In this case, 
the shocked gas adiabatically expands and cools after the 
shock wave has traversed the dense wind, and the gas is 
ultimately expected to return to free expansion. We do 
not examine the details of this process here, but note 
that energy conserving interaction with a wind appears 
not to be ruled out for the Crab. 


3.3. Detailed Pulsar Nebula Model 

In the previous sections, we assumed a simple form for 
the freely expanding supernova density profile as well as 
steady pulsar power input to enable analytic estimates of 
the model parameters. Here, we use a s upernova density 
profil e based on detailed calculations (Matzner & McKee 
1999) and the standard evolution of the pulsar power 
(equation [3]). The ordinary differential equations 
(18) for the evolution of the pulsar bubble are solved with 
aKunga-Kutta scheme. 

In previous discussions of pulsar nebula expansion in a 
supernova, the supernova has a normal energy (10 51 ergs) 
and the expansion of the nebula is just in the inner flat 
part of the density profile. For a low energy supernova, 
the transition to the outer steeply dropping part of the 
profile may become important. The supernova density 
profile is shown in Figure [2] for an energy of 1 x 10 50 ergs 
and an ejecta mass of 5 M©. The model is intended for 
a red supergiant prog enitor star, as might be expected 
(iTomi naga et al.|2013). The model of |Matzner & McKee 
(1999J lor red supergiants has a parameter g, which is 
related to the core-envelope mass ratio and is in the range 
0.3 — 0.8; we choose q = 0.5. At small velocities, the 
profile approaches a flat power law with index S = 0.21 
and at high velocities a steep power law with index p = 
11.7 (Figure [2|. 

A problem with comparing the model to the actual 
Crab Nebula is that the model assumes that the pul¬ 
sar bubble sweeps up a thin shell, while the Crab fil¬ 
aments are spread in radius, presumably as a result of 
the Rayleigh-Taylor instability. We consider two cases 
to cover the possibilities. In the first, the diameter of 
the Crab is taken to be 3.3 pc, consistent with the outer 
boundary of the Crab. Assuming that the filaments at 
the outer boundary have the typical At = 75 yr, we have 
a velocity of 1840 km s -1 . For the other case, we assume 
that the velocity of the dense shell is characterized bv 
that of the dense filaments, about 1260 km s -1 (Temim 


et al. 2006). Now assuming that At = 75 yr leads to a 
diameter of 2.23 pc. We expect that the two cases will 
bracket the actual Crab Nebula. 

In order to estimate the goodness of fit of assumed 
models, we calculated 


x 2 = [(Rfit ~ Ro)/Ro} 2 + [( v fit - Vo)/Vo] 2 (18) 

where Rfit and Vfu are the estimated observed val¬ 
ues of the shell radius and velocity, and Ro and Vo 
are the results of models for various values of ejecta 
mass Msn and supernova energy Es n- In Figures [3] 
and 01 the darkness of the plots indicates the value of 
log(V). Figure I shows models with Rfit — 1-65 pc 
and Vfu = 1840 km s _1 , where the model parame¬ 
ters are the total ejecta mass, Msn , and the super¬ 
nova energy Esn • If Msn is 6.6 i 1 M©, or less, as 



0.5 1.0 1.5 2.0 


Esn (erg) xlO 50 

Fig. 3. — Results of models as a function of explosion energy 
and total ejecta mass. The shading gives the chi-squared value in 
fitting a shell radius of 1.65 pc and velocity of 1840 km s -1 . 
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0 

Fig. 4. — Same as Figure 2, but fitting a shell radius of 1.1 pc 
and velocity of 1260 km s -1 . 
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discussed in Section 2, the energy of the supernova is 
Esn ~ 1.0 zb 0.5 x 10 50 ergs. Figure [4] shows results for 
the assumption of the compact shell with Rf it = 1.1 pc 
and Vfu = 1260 km s -1 . The acceptable models occur 
at low energy with Esn < 0.5 x 10 50 ergs. 

Another constraint comes from the swept up mass in 
the shell, whic h is estim ated as 4.6 ±1.8 M© from the ob¬ 
servations (Fesen et al.j[l997|. Figure [ 5 ] shows the swept 
up mass resulting from a given total ejecta mass, Msn , 
and the supernova energy Esn- With the limit on Msn , 
a low value for the energy, < 10 50 ergs, is indicated, 
in agreement with the approximate model described in 
Section 2.1. A recent esti mate of the mass in th e Crab 
filaments is 7.2 d= 0.5 M© (|Owen & Barlow||2015|); using 
this higher mass would yield an even stricter limit on the 
energy. 

Also, the density profile of the freely expanding 
gas near the edge of the Crab Nebula (velocity ~ 
1600 km s -1 ) has a density gradient between that of the 
outer steep power law and the inner fla t region (Figure 


This is in accord with the findings of 


Sollerman et al. 


(zOOO) on the distribution of C IV outside the main neT> 
ula. A higher energy supernova would give a relatively 
flat density profile for the supernova ejecta just outside 
the main body of the Crab. 


4. CONCLUSIONS 
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Fig. 5. — The swept up shell mass resulting from a given explosion 
energy and total ejecta mass. The labels on the solid lines give the 
swept up mass in M©. The black li nes show 4 .6 d= 1.8M© where 
1.8 M© is one standard deviation (|Fesen et al.j|[1997|); the blue lines 
show +0.5, —0.5, —1.5, —2.0 standard deviations, respectively. 

Independent lines of evidence support a low energy 
10 50 erg) for the supernova giving rise to the Crab 


Nebula. The evidence includes the radius and accelera¬ 
tion of the Crab shell, the mass of swept up ejecta in the 
Crab shell, the presence of a radiative shock wave around 
part of the Crab shell, and the decreasing amount of C 
IV in the region surrounding the nebula. 

The model proposed here combines the low energy su¬ 
pernova view of the Crab supernova with expansion into 
freely expanding gas. However, we do not consider the 
interaction of the freely expanding gas with its surround¬ 
ings. The nature of the presupernova mass loss and the 
interaction of the supernova with the mass loss remains 
unclear. Further studies of the region around the Crab 
Nebula are warranted. 
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